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O . Abstract 

^ The numerical solutions of nonlocal and local Boltzmann kinetic equations for the 

simulation of peripheral and central heavy ion reactions are compared. The experimental 
finding of enhancement of mid-rapidity matter shows the necessity to include nonlocal 
I corrections. While the in-medium cross section changes the number of collisions and not 

the transferred energy, the nonlocal scenario changes the energy transferred during colli- 
sions. The renormalisation of quasiparticle energies by using the Pauli-rejected collisions 
ly^ I leads to a further enhancement of mid-rapidity matter distribution and is accompanied 

' by a dynamical softening of the equation of state. 

The simulation results are parameterised in terms of time dependent thermodynamical 
I variables in the Fermi liquid sense. This allows one to discuss dynamical trajectories in 

phase space. A combination of thermodynamical observables is constructed which locates 
instabilities and points of possible phase transition under iso-nothing conditions. Two 
different mechanisms of instability, a short time surface-dominated instability and later a 
spinodal-dominated volume instability is found. The latter one occurs only if the incident 
energies do not exceed significantly the Fermi energy and might be attributed to spinodal 
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, decomposition. 
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^ : 1 Introduction 

The formation of a neck-like structure in peripheral heavy ion reactions and the impact on 
the fragmentation mechanism and production of light charged particles has been discussed for 
a couple of years [Q- Theoretical investigations suggest that the neck is not formed in 
usual heavy ion simulations starting from the Landau equation ^, H or BUU equations H, 
including additional mean field fluctuations derived in 0, ^ and tested |10]. The inclusion of 



fluctuations in the Boltzmann (BUU) equation has been investigated resulting in Boltzmann- 
Langevin pictures |]1T], |12|, 



We will take here the point of view that the fluctuations should arise by themselves in a 
proper kinetic description where all relevant correlations are included in the collision integral. 
The collision will then cause both a dephasing and fluctuation by itself. This procedure without 
additional assumption about fluctuations has been given by the nonlocal kinetic theory [l^, |16 



171 and applied to heavy ion collisions in [Q- [0. We claim that the derived nonlocal off-set 
in the collision procedure induces fluctuations in the density and consequently in the mean-field 
which are similar to the one assumed in Boltzmann Langevin approaches above. 

Recent INDRA observation shows an enhancement of emitted matter in the region of almost 
zero relative velocity which means that matter is stopped during the reaction and stays almost 
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at rest |21], ^ Q . This enhancement of mid-rapidity distribution can possibly be associated 
with a pronounced neck formation of matter. 

We use for the description the nonlocal kinetic equation ^ for the one - particle 
distribution function 
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with Enskog-type shifts of arguments |5|, |2|]: /i = /a(k,r,t), /a = /^(p, r - A2, t), /a = 
/a(k — q— Ax, r — A3, t — At) , and = /;,(p + q— A^^-, r — A4, t — A^) . The effective scattering 
measure, the T-matrix is centred in all shifts. The quasiparticle energy e contains the mean 
field as well as the correlated self energy. In agreement with p5| , all gradient corrections 
are given by derivatives of the scattering phase shift = Im InT^f k, p, q,t, r). 
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After derivatives, A's are evaluated at the energy shell Q £3 + 6:4. Neglecting these shifts the 
usual BUU scenario appears. The scattering integral of the non-local kinetic equation derived 
in 1^ ^ corresponds to the picture of a collision as seen in Figure 

Despite its complicated form it is possible to solve this kinetic equation with standard Boltz- 

27| , p!9| . To this end we have calculated 
These shifts and the modifications of 
standard BUU or QMD code are available from the author. The numerical solution of the 



mann numerical codes and to implement the shifts 
the shifts for different realistic nuclear potentials flS 



nonlocal kinetic equation has shown an observable effect in the dynamical particle spectra |16 
as well as in the charge density distribution [|19|. The high energetic tails of the spectrum 
are enhanced due to more energetic two-particle collisions in the early phase of nuclear colli- 
sion. Therefore the nonlocal corrections lead to an enhanced production of preequilibrium high 
energetic particles. 

Besides the nonlocal shifts and cross section which has been calculated from realistic poten- 
tials, the interaction affects, however, the free motion of particles between individual collisions. 
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The dominant effect is due to mean-field forces which bind the nucleus together, accelerat- 
ing particles close to the surface towards the centre. These mean-field forces are conveniently 
included via potentials of Skyrme type 

e = f^A(^)^B(^Y. (3) 

For a derivation of collision integrals and the Skyrme potential from the same microscopic 
footing, see |p8 |. 



Beside forces, the interaction also modifies the velocity with which a particle of a given 
momentum propagates in the system. This effect is known as the mass renormalisation. A 
numerical implementation of the renormalised mass is rather involved since a plain use of the 
renormalised mass instead of the free one leads to incorrect currents. Within the Landau 
concept of quasiparticles, this problem is cured by the back flow, but it is not obvious how to 
implement the back flow within the BUU simulation scheme. In our studies, we circumvent the 
problem of back flows using explicit zero-angle collisions to which we add a non-local correction. 
One can show that indeed the non-local shifts create just the dynamical mass renormalisation 
if used for the Pauli-rejected events. The incorporation of a nonlocal jump without performing 
finite angle collisions for such events correspond exactly to the quasiparticle and dynamical 



mass renormalisation 19 



Finally, we would like to comment on properties of the proposed simulation scheme. The 
renormalisation depends on the distribution of particles in surrounding medium. It has four 
nice properties: (i) the renormalisation vanishes as the local density goes to zero, (ii) the 
renormalisation vanishes when a high temperature closes the Luttinger gap because all collisions 
will be at finite angles, (iii) the anisotropy of the quasiparticle velocity in a presence of a non- 
zero current in medium is automatically covered, and (iv) the backfiows connected to the mass 
renormalisation are covered because both particles jump keeping the centre of mass fixed. Last 
but not least, the simulation does not require to introduce new time-demanding procedures, 
one can simply use the scattering events which are merely rejected in standard simulation codes 
by Pauli-blocking. 



2 Numerical results 

Let us discuss the proposed correction to the local and ideal (no quasiparticle renormalisation) 
Boltzmann (BUU) simulation. First we introduce the pure nonlocal corrections and then we 
discuss the quasiparticle renormalisation. 

The evolution of the density can be seen in the corresponding left pictures of figure H for the 
BUU (left panel) and nonlocal scenario (middle panel) as well as the additional quasiparticle 
renormalisation (right panel). We see that the nonlocal scenario leads to a longer and more 
pronounced neck formation between 200 — 240fm/c while the BUU breaks apart already at 
200fm/c. 

The question arises whether this pronounced neck formation and increase of mid-rapidity 
distribution is simply occurring by more collisions and corresponding correlations as called in- 



medium effect. This is however only the case for smaller impact parameters [^]. We see in 
the next figure ^ the number of collisions per time for the simulation where in the local BUU 
scenario the cross section has been doubled. The number of collisions is visibly enhanced by 
doubling the cross section while for the nonlocal scenario we get only a slight enhancement at 
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Figure 2: First part. 
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Figure 2: The time evolution of jf^Ta +7^^ Au collisions at Eiab/A = 33 MeV and 8 fm impact 
parameter in the local BUU (left), non-local BUU (middle), and the non-local BUU with 
quasiparticle renormalisation (right). Plots in the first column show the {x — 2;)-density cut 
of 30 fmx30 fm where Ta as projectile comes from below. The mass momenta are shown by 
arrows. The corresponding second columns give the charge density distribution versus relative 
velocity in cm/ns where the target like distribution of Au is on the left and the projectile like 
distributions of Ta on the right. 
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the beginning and later even lower values with respect to local BUU. The latter fact comes 
from the earlier decomposition of matter in the nonlocal scenario. 
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Figure 3: The time evolution of the num- 
ber of nucleon collisions for Ta + Au at 
Eiab/A = 33 MeV and different impact 
parameter in the BUU (thick black line), 
non-local kinetic model (broken line) and 
the local BUU with a doubled cross sec- 
tion (thin dark line). The impact param- 
eter is 8 fm. 



The corresponding transverse and kinetic energies in figure ^ (8fm impact parameter) show 
that the transverse and longitudinal energy is almost not changed compared with local BUU. 
Oppositely the nonlocal scenario leads to an increase of transverse energy of about 2MeV and 
about IMeV in longitudinal energy. We conclude that the increase of cross section leads to 
a higher number of collisions but not to more dissipated energy while the nonlocal scenario 
does not change the number of collisions much but the energy dissipated during the collisions. 
Returning to the discussion of pronounced neck formation in figure ^ above we see now that 
the quality rather than the quantity of collisions is what produces the neck. 
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Figure 4: The time evolution of the longi- 
tudinal (thin lines) and transverse energy 
(thick lines) including Fermi motion of nu- 
cleon collisions for Ta+Au at Eiab/A = 33 
MeV in the BUU (black line), nonlocal 
kinetic equation (light line) and the lo- 
cal BUU but with twice the cross section 
(medium line). 



Now we can proceed and discuss the charge matter distribution with respect to the velocity. 
We plotted in figure | also the normalised charge distribution versus velocity and see that after 
160fm/c we have an appreciable higher mid-rapidity distribution for the nonlocal scenario (mid 
panel) than the BUU (left panel). Together with the observation that for nonlocal scenario we 
have a pronounced neck formation we see indeed that the neck formation is accompanied with 
high mid- velocity distribution of matter. We see in figure |^ (right panel) that the mid-rapidity 
distribution of matter is once more enhanced for quasiparticle renormalisation in comparison 
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to the nonlocal scenario without quasiparticle renormalisation. The detailed comparison of the 
time evolutions of the transverse energy for 8fm impact parameter can be seen in figure ^ 
We recognise that the transverse energies including quasiparticle renormalisation are similar to 
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Figure 5: The time evolution 
of the transverse energy in- 
cluding Fermi motion for Ta + 
Au at Eiab/A = 33 MeV 
and 8fm impact parameter in 
the BUU (black line), nonlo- 
cal kinetic equation (dashed 
line), the local BUU but twice 
cross section (dashed dotted 
line) and the nonlocal scenario 
with quasiparticle renormalisa- 
tion (long dashed line). 



the nonlocal scenario and higher than the BUU or BUU with twice the cross section. However 
please remark that the period of oscillation in the transverse energy which corresponds to a giant 
resonance becomes larger for the case with quasiparticle renormalisation. Since therefore the 
energy of this resonance decreases we can conclude that the (nuclear) compressibility has been 
decreased by the quasiparticle renormalisation. Sometimes this quasiparticle renormalisation 
has been introduced by momentum dependent mean-fields. The effect is known to soften 
the equation of state. We see here that we get a dynamical quasiparticle renormalisation 
and a softening of equation of state. This softening of equation of state is already slightly 
remarkable when the nonlocal scenario is compared with BUU. With additional quasiparticle 
renormalisation we see that this is much pronounced. 

We want to repeat that the dynamical quasiparticle renormalisation which leads to a soft- 
ening of the equation of state enhances the mid rapidity distribution. In contrast a mere soft 
static parametrisation of the mean-field does not change the mid-rapidity emission appreciably 



2.1 Comparison with experiments 

The BUU simulations will now be compared to one experiment performed with INDRA at 
GANIL, the Ta + Au collision at Eiab/A = 33 MeV. For the identification with experimental 
data we select events which show a clear one fragment structure. This correspond to events 
where we have clear target and projectile like residues. Since the used kinetic theory is not 
capable to describe dynamical fragment formation we believe that these events are the ones 
which are at least describable within our frame. Next we use impact parameter cuts with 
respect to the transverse energy since this shows in all simulations a fairly good correlation 
to the impact parameter. In our numerical results we see almost linear correlations between 
impact parameter, maximal velocity and the convenient ratio between transverse and total 
kinetic energy as seen in figure ^. 
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Figure 6: The experimental 
charge distribution of matter 
(dotted hne) versus velocity in 
comparison within the BUU 
(thin solid line) and the non- 
local model with quasiparticle 
renormalisation (thick line) . 
The maximum velocity versus 
ratio of longitudinal to total ki- 
netic energy of Ta + Au colli- 
sions at Eiab/A = 33 MeV is 
given below. The selected ex- 
perimental cuts are given by 
dots. 
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For each selected experimental transverse energy bin we can plot now the maximum veloc- 
ity versus the ratio of the transverse to kinetic energy. We see in figure |^ that the numerical 
velocity damping agrees with the experimental selection only for very peripheral collisions. For 
such events we plot in the figure ^ the charge density distribution and compare the experiment 
with the simulation. These charge density distributions have been obtained using the procedure 
described in reference The Data are represented by light grey points, the standard BUU 
calculation by the thin line and the non-local BUU with quasiparticle renormalisation calcu- 
lation by the thick line. A reasonable agreement is found for the nonlocal scenario including 
quasiparticle renormalisation while simple BUU fails to reproduce mid-rapidity matter. 



3 Nonequilibrium Thermodynamics 

In order to understand the different thermodynamics induced by nonlocal corrections and there- 
fore virial correlations, we multiply the kinetic equation with l,p,e and obtain the balance for 
the particle density n, the momentum density J and the energy density S, see details in 



Please note that besides the mean field (^ we have also a Born correlation term, see |2^, |2S 

The correlational parts of the density, pressure and energy are arising from genuine two- 
particle correlations beyond Born approximation which are also derived from the balance equa- 
tions of nonlocal kinetic equations |T5|, E3. It has been shown that they establish the complete 



conservation laws. While these correlated parts are present in the numerical results and can be 
shown to contribute to the conservation laws we will only discuss the thermodynamical prop- 
erties in terms of quasiparticle quantities to compare as close as possible with the mean field or 
local BUU expressions. The discussion of these correlated two - particle quantities is devoted 
to a separate consideration. 

From the distribution function /(p,r, t) the local density, current and energy densities are 



8 



given by 

/ n(r,t) 
J(r,t) 



dp ^ ^ 



P /(P,r,t) (4) 



m 



which are computed directly from the numerical solution of the kinetic equation in terms of 
test particles. Please note that the above kinetic energy includes the Fermi motion. 

The global variables per particle number like kinetic energy, Fermi energy and collective 
energy are obtained by spatial integration 

£f[nir,t)]n{r,t) / / drn{r,t) (5) 
J(r,t)Vmn(r,t)/ 



where we have used the local density approximation |^ . Now we adopt the picture of Fermi 
liquid theory which connects the temperature with the kinetic energy as 

SKit) = ^Sp{t) + S^t) + ^^fitf (6) 

from which we deduce the global temperature. The definition of temperature is by no means 
obvious since it is in principle an equilibrium quantity. One has several possibilities to define 
a time dependent equivalent temperature which should approach the equilibrium value when 
the system approaches equihbrium. In the definition of slope temperatures has been 

discussed and compared to local space dependent temperature fits of the distribution function 
of matter. This seems to be a good measure for higher energetic collisions in the relativistic 
regime. Since we restrict us here to collisions in the Fermi energy domain and do not want to 
add coalescence models we will not use the slope temperature. Moreover we define the global 
temperature in terms of global energies which are obtained by local quantities rather than by 
defining a local temperature itself. This has the advantage that we do not consider local energy 
fluctuations but only a mean evolution of temperature. 
The mean field part of the energy is given by 

u(t) = £.pw-£.W-£-W = /<<.(^^ + b;^)//*«(M) (7) 

from which one deduces the pressure per particle 

Pit) = liSKit) - £Mt)) + + ldr (^A^ + i?^^^)// drnir,t). (8) 

In order to compare now the local BUU with the nonlocal BUU scenario we consider the 
energy which would be the total energy in the local BUU without Coulomb energy 

S{t)=SK{t)~S,oii{t) + U{t). (9) 

This expression does not contain the two - particle correlation energy which is zero for BUU 
and the Coulomb energy. The reason for considering this energy for dynamical trajectories is 
that we want to follow the trajectories in the picture of mean field and usual spinodal plots. 

To define the density exhibits to some extend a problem. Depending on the considered 
volume sphere we obtain different global densities. We follow here the point of view that the 
mean square radius will be used as a sphere to define the global density. This is also supported 
by the observation that the mean square radius follows the visible compression. 
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3.1 Iso-nothing conditions in equilibrium 

Let us first recall the figures of mean field isotherms in equilibrium. We obtain the typical 
van der Waals curves in figure 0. Since we have neither isothermal nor isochoric nor isobaric 
conditions in simulations, shortly since we have iso-nothing conditions, we have to find a rep- 
resentation of the phase transition curves which are independent of temperature but which 
reflect the main features of phase transitions. This can be achieved by the product of energy 
and pressure density versus energy density in figure |^ below. This plot shows that all instable 
isotherms exhibit a minimum in the left lower quarter. There the energy is negative denoting 
bound state conditions but the pressure is already positive which means the system is unstable. 
The first isotherm above the critical one does not touch this quarter but remains in the right 
upper quarter where the energy and pressure are both positive and the system is expanding 
and decomposing unboundly. The left upper quarter denotes negative pressure and energy 
indicating that the system is bound and stable. 

In order to achieve now a temperature independent plot we scale both axes of figure |^ 
(below) with a temperature dependent polynomial and achieve that all critical isotherms are 
collapsing on one curve in the left lower quarter. The first isotherms above the critical one does 
not enter the left lower quarter. We consider this scaling as adequate for iso-nothing conditions. 
A phase transition should be possible to observe if there occurs a minimum in the left half of 
this plot at negative energies. 





Figure 7: The isotherms for the pressure 
density versus volume and for the prod- 
uct of pressure and energy density ver- 
sus energy density (above). The temper- 
atures are T = 1, 4, 7, 10, 13, 16, 19MeV. 
The isotherms of the product of pressure 
and energy density versus energy density 
scaled by a temperature polynomial f(T), 
g{T) (below). All critical isotherms col- 
lapse on one line in the left lower quarter. 



The idea of plotting combinations of pressure and energy is similar to the one of softest 
point [^] in analysing QCD phase transitions. There the simple pressure over energy ratio 
leads to a temperature independent plot due to ultrarelativistic energy - temperature relations. 
In our case we have a Fermi liquid behaviour at low temperatures and the scaling temperature 
dependent polynomials, / and g, can be found in ||2^ which are producing a temperature 
independent plot in figure |^. 
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Figure 8: The dynamical trajec- 
tories of the energy (1), density 
and temperature in the nonlocal 
(gray thick) and in the local BUU 
(black thin) scenario. The consid- 
ered reaction is ^^^Xe on ^^^S'n at 
25MeV lab energy. The dots mark 
the times in steps of 20fm/ c up to to- 
tal of 300fm/c. To guide the eye the 
zero temperature mean field energy 
(thick line) and the pressure (thin 
line) is plotted in the upper left pic- 
ture and in the right figures the spin- 
odal line for infinite matter is given. 
The scaled combinational plot anal- 
ogous to figure ^ is given in the left 
lower plane. 

3.2 Simulation results 

Let us now inspect the dynamical trajectories for the above defined temperature, density and 
energy. In figure ^ we have plotted the dynamical trajectories for a charge - symmetric reaction 
of Xe on Sn at 25MeV lab - energy. The solution of the nonlocal kinetic equation is compared 
to the local BUU one. One sees in the temperature versus density plane that the point of 
highest compression is reached around 60fm/c with a temperature of 9 MeV. 

After this point of highest overlap or fusion phase we have an expansion phase where the 
density and temperature is decreasing. While the compression phase is developing similar for 
the BUU and for the nonlocal kinetic equation we see now differences in the development. First 
the temperature of the nonlocal kinetic equation is around 2MeV higher than the local BUU 
result. This is due to the release of correlation energy into kinetic energy which is not present 
in the local BUU scenario. After this expansion stage until times of 120fm/c we see that the 
BUU trajectories come to a rest inside the spinodal region while the nonlocal scenario leads to 
a further decay. This can be seen by the continuous decrease of density and increase of energy. 
Since matter is more decomposed with the nonlocal kinetic equation we also heat the system 
more due to Coulomb acceleration. This leads to the enhancement of temperature compared to 
BUU. An oscillating behaviour occurs at later times which reflects an interplay between short 
range correlation and long range Coulomb repulsion. The decomposition leads almost to free 
gaseous matter after 300fm/c as can be seen in the energy versus density plot. 

Please note that although the trajectories seem to equilibrate inside the spinodal region 
when one considers the temperature versus density plane, we see that in the corresponding 
energy versus temperature plane the trajectories travel already outside the spinodal region. 
This underlines the importance to investigate the region of spinodal decomposition in terms 
of a three dimensional plot instead of a two dimensional one like in the recently discussed 
caloric curve plots. Different experimental situations lead to different curves as long as the 
third coordinate (pressure or density) remains undetermined. 

The iso-nothing plot analogous to the left lower plot of figure shows that the point of 
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highest compression is hnked to a first instabihty seen as a pronounced minimum of the trajec- 
tory in the left quarter. If we do not exceed the Fermi energy domain, this is connected with a 



pronounced surface emission and connected with anomalous velocity profiles |3J]. We will call 
this phase surface emission instability further on. At 180 fm/c we see a second minimum which 
is taking place inside the spinodal region. This instability we might now attribute to spinodal 
decomposition since the trajectories developing slower and remain inside the spinodal region. 
The BUU shows the same qualitative minima but the matter rebounds and the trajectories 
move towards negative energies again. In opposition the nonlocal scenario leads to a further 
decomposition of matter as described above. 

In the next figure ^ we have plotted the same reaction as in figure |^ but at higher energy 
of 33MeV. We recognise a higher compression density and temperature than compared to the 
lower bombarding energy. Consequently the trajectories develop further towards the unbound 
region of positive energy after 300fm/c. While the first surface emission instability is strongly 
pronounced we see that the second minimum in the iso-nothing plot is already weaker indicat- 
ing that the role of spinodal decomposition is diminished. The trajectories in the temperature 
versus density plot come still in the spinodal region at rest but travel already outside the spin- 
odal region if the energy versus temperature plot is considered. This shows that the trajectories 
start to develop too fast to suffer much spinodal decomposition. Oppositely, at energies around 
the Fermi energy spinodal decomposition might be possible and is presumably connected with 



anomalous velocity and density profiles |^4[. If we now plot the same reaction at 50MeV in 
figure |lO| we see that the trajectories come at rest outside the spinodal region whatever plot 
is used and no second minima is seen anymore in the iso-nothing plot. But, the surface emis- 
sion instability is still very pronounced and is probably here the leading mechanism of matter 
disintegration. We might now search for a situation where we have the opposite extreme that 
is we search for a reaction with as less as possible surface emission instability and as much 
as possible spinodal decomposition. For this reason we might think on asymmetric reactions 
since the different sizes of the colliding nuclei might suppress the surface emission. Indeed as 
can be seen in figure |TT] for an asymmetric reaction of Ni on Au at 25MeV lab-energy with 
nearly the same total charge as in the reaction before that the surface emission instability is 
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less pronounced while the spinodal instability is much more important. There appears even a 
third minima showing that the matter suffers spinodal decomposition perhaps more than once 
if the bombarding energy is low enough and a long oscillating piece of matter is developing. 

The higher bombarding energies now show the same qualitative effect which pronounces the 
surface emission instability and reduces the importance of the spinodal decomposition. Much 
smaller compression densities and temperatures are reached in these reactions compared to the 
more symmetric case of Xe on Sn. 



4 Summary 

The extension of BUU simulations by nonlocal shifts and quasiparticle renormalisation has been 
presented and compared to recent experimental data on mid rapidity charge distributions. It 
is found that both the nonlocal shifts as well as the quasiparticle renormalisation must be 
included in order to get the observed mid-rapidity matter enhancement. 

The inclusion of quasiparticle renormalisation has been performed by using the normally 
excluded events by Pauli blocking. Since the quasiparticle renormalisation and corresponding 
effective mass features can be considered as zero angle collisions they can be realized by nonlocal 
shifts for the scattering events which are normally rejected. This means that one has to perform 
the advection step for the cases of Pauli blocked collisions without colliding the particles. 
Besides giving a better description of experiments, this has the effect of a dynamically softening 
of equation of state seen in longer oscillations of giant compressional resonance. 

In this way we present a combined picture including nonlocal off-sets representing the non- 
local character of scattering, which leads to virial correlations with the quasiparticle renormal- 
isation, and as a result to mean field fluctuations. We propose that no additional stochasticity 
need to be assumed in order to get realistic fluctuations. 

The nonlocal kinetic theory leads to a different nonequilibrium thermodynamics compared 
to the local BUU. We see basically a higher energetic particle spectra and a higher transversal 
temperature of 2MeV. This is attributed to the conversion of two-particle correlation energy 
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Figure 11: The same dynamical tra- 
jectories as in figure but for a re- 
action ^^Ni on i^Mm at 25MeV lab 
energy. 
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Table 1: The prediction of the leading mechanisms of matter disintegration for two reactions 
with equal total charge but symmetric and asymmetric entrance channels. Surface compression 
is denoted by C and spinodal decomposition by S. 



into kinetic energy which is of course absent in local BUU scenario. 

By constructing a temperature independent combination of thermodynamical variables we 
are able to investigate the signals of phase transitions under iso-nothing conditions. Two 
mechanisms of instability have been identified: surface emission instability and spinodal de- 
composition. We predict for the currently investigated reactions seen in table |I| which effect 
should be the leading one for matter decomposition. In the reactions with bombarding energies 
higher than the Fermi energy the fast surface eruption happens outside the spinodal region. 
For even higher energies there is not enough time for the system to rest at the spinodal region. 
The trajectories simply move through the spinodal and the system decays before it comes to 
an equilibrium - like state inside the spinodal region. Around the Fermi energy the spinodal 
decomposition might occur and is accompanied by an anomalous velocity and density profile 
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